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Abstract 

Helix-coil transitions in poly-alanine molecules of length 10 are studied by 
multicanonical Monte Carlo simulations. The solvation effects are included 
by either a distance-dependent dielectric permittivity or by a term that is 
proportional to the solvent-accessible surface area of the peptide. We found a 
strong dependence of the characteristics of the helix-coil transition from the 
details of the solvation model. 
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1 Introduction 

There has been recently a renewed interest in the conditions under which a-helices, 
a common structure in proteins, are formed or dissolved. It is long known that a- 
helices undergo a sharp transition towards a random coil state when the temperature 
is increased. The characteristics of this so-called helix-coil transition have been 
studied extensively most recently in Refs. [Q, [|. In Refs. ||, |5| evidence was 
presented that the helix-coil transition in polyalanine exhibits a true thermodynamic 
phase transition when interactions between all atoms in the molecule are taken into 
account @, ^ . 

The later results were obtained from gas-phase simulations of poly-alanine. While 
there is some experimental evidence || supporting the numerical results of these gas- 
phase simulations, the question remains how these results relate to the biologically 
more relevant case of solvated molecules. First investigations of this question were 
described in Refs.R, ||] where it was claimed that the transition temperature is lower 
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in water than in vacuum. However, that investigation relies on a single represen- 
tation of the protein-water interaction and the dependence of their results on the 
details of the solvation term is not clear. 

In this paper, we have investigated how the characteristics of helix-coil transition 
change with the details of the solvation term. For this purpose, we have performed 
multicanonical simulations of polyalanine molecules of length 10. The protein-water 
interaction was included in two ways: either by a distance-dependent dielectric per- 
mittivity or by a term that is proportional to the solvent- accessible surface area of 
the peptide. For the later case we have considered four different parameter sets: 
OONS @, JRF |3, W92 Q and SCH pf. Quantities such as the energy, helicity 
and susceptibility were calculated as function of temperature. Our result were com- 
pared with that of gas phase simulations. A strong dependence of the characteristics 
of the helix-coil transition from the details of the solvation term was found. 



2 Methods 



Our investigation of the helix-coil transition for polyalanine is based on a detailed, 
all-atom representation of that homopolymer. The interaction between the atoms 
was described by a standard force field, ECEPP/2, ]nj (as implemented in the pro- 



gram package SMMP [14]]) and is given by: 
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Here, (in A) is the distance between the atoms i and j, and xi is t ne ^"th torsion 
angle. We have chosen ECEPP/2 instead of the newer ECEPP/3 differs because 
this choice allows a more easy comparison with our previous work. Both force fields 
differ from each other only in the way in which prolines and end groups are treated. 
In preliminary polyalanine simulations we found no qualitative differences in our 
results when ECEPP/3 was used instead of ECEPP/2 (data not shown). 
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The interactions between our homo-oligomer and water are approximated by 
means of two implicit water models. In the first model (DDE) the electrostatic 
interactions in the presence of water rely on a distance dependent electrostatic per- 
mittivity [15]: 



e(r) =D- [(srf + 2sr + 2}e~ sr . (6) 

For the parameters D and s empirical values are chosen such that for large distances 
the permittivity takes the value of bulk water (e fh 80), and the value e = 2 for short 
distances (protein interior space). Equation (^) is the result of interpolation of two 
types of interactions. For short distances it models the interaction of two charges 
placed in continuum medium, while over long distances it represents a Debye curve. 
This is clearly a gross over-simplification of protein-solvent interactions. However, 
approximating solvation effects by a distance-dependent dielectric permittivity was 



used by many authors to study the proteins and nucleic acids (e.g. ||16|| ) since it 
does not significantly slow down protein simulations below that of simple in vacuo 
simulations. 

In another common approximation of the protein-solvent interaction one assumes 
that the free energy contributions from atomic groups immersed in the protein in- 
terior differ from contributions of groups exposed to the water. It is commonly 



accepted || [11], [L7|, [18[] that this free-energy difference is proportional to the surface 
area of the atomic group which is exposed to the solvent. Within this approxima- 
tion, the total solvation energy of a protein is given by the sum of contributions 
from each solvated atomic groups: 

E sot = J2o-iAi, (7) 

i 

where E so i is the solvation energy, A t is the conformational dependent solvent acces- 
sible area of the surface of the i — th atom and <7j is the atomic solvation parameter 
for the atom i. The summation is extended over all atomic groups. The solvation 
parameters are evaluated experimentally by measuring the free energy needed to 
bring the group from a nonpolar environment (usually octanol or ethanol are used 
as convenient compounds) into water. Many sets of solvation parameters were eval- 
uated by several authors with different methods, but unfortunately it is not always 
obvious which one is the most appropriate one. The sets we study here are named by 
us OONS §, JRF 0, W92 fTJ and SCH [TJ, and are described in the respective 
references. 
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Simulations of detailed models of biological macromolecules are notoriously dif- 
ficult. This is because the various competing interactions within the polymer lead 
to an energy landscape characterized by a multitude of local minima. Hence, in 
the low-temperature region, canonical Monte Carlo or molecular dynamics simu- 
lations will tend to get trapped in one of these minima and the simulation will 
not thermalize within the available CPU time. Only recently, with the introduc- 
tion of new and sophisticated algorithms such as multicanonical sampling |19| and 
other generalized- ensemble techniques [p(J was it possible to alleviate this problem 
in protein simulations [^TJ. For polyalanine, both the failure of standard Monte 



Carlo techniques and the superior performance of the multicanonical algorithm are 



extensively documented in earlier work [22]. For this reason, we use again this 



sophisticated simulation technique for our project. 

In the multicanonical algorithm [|19| conformations with energy E are assigned 
a weight w mu (E) oc l/n(E). Here, n(E) is the density of states. A simulation with 
this weight will lead to a uniform distribution of energy: 

P mu {E) oc n(E) w mu (E) = const . (8) 

This is because the simulation generates a ID random walk in the energy space, 
allowing itself to escape from any local minimum. Since a large range of energies are 
sampled, one can use the reweighting techniques [^3| to calculate thermodynamic 
quantities over a wide range of temperatures T by 

Jdx A{x) w-\E{x)) e- pE W 
< T Jdx w-i(E(x)) e-M*) ' ( } 

where x stands for configurations. 

Unlike in the case of canonical simulations the weights 

W (E) = n-\E) = e~ s{E) (10) 

are not a priori known. Instead estimators for these weights have to be determined. 
This is often done by an iterative procedure in which for reasons of numerical sta- 
bility Eq. [H] is replaced by 



The multicanonical parameters /3(E) and a{E) are defined through 

a/JP . S{E')-S{E) f , E : 

m = E'-E and a{E) = I a(E') + ((3(E>) - (3(E))E' ,E<E„ 



max 
max 



(12) 
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with E and E' adjacent bins in the array S(E). The (3(E) are then iteratively 
updated by the relation 



(3 l+l (E) = (3\E) + g (E) ■ (\nH\E>) - \nH\E)) /{Ef - E) , (13) 

in which H l (E) is the histogram of the i-th run (and H(E) > 1). In Ref. |24| the 
factor go(E) in Eq. [13] was defined through 



a ° {E)= ¥mr^hm with g,(£) = Iw+ffl • (14) 



The above relation assumes that the histogram H(E) counts independent events 
which is in general not true. Hence, it is more appropriate and leads to a faster 
convergence of /3(E) if the array g l (E) in Eq. ^ is instead defined by 

y v ; K*{E') + K*(E) v ; 

where the auxiliary array K(E) now counts only the number of independent visits 
at energy E. 

With the above described iterative procedure we needed 200,000 sweeps for the 
weight factor calculations. All thermodynamic quantities were then estimated from 
one production run of 1, 000, 000 Monte Carlo sweeps starting from a random initial 
conformation, i.e. without introducing any bias. 

3 Results and Discussion 

In previous gas-phase simulations of poly-alanine 0, ||], [5], |22| we observed at T = 
430 K a pronounced transition between a high-temperature phase dominated by dis- 
ordered coil structures and an ordered phase with single, extended helices. A natural 
order parameter for this helix-coil transition is the average number < nu(T) > of 



residues in the oligomer which are part of an a— helix. Following earlier work \T1 
we define a residue as helical if the pair of backbone dihedral angles 0, ip takes a 
value in the range (— 70±20, — 37±20). In Fig. la this order parameter is displayed 
as function of temperature for a gas-phase simulation (GP) of Alaio and simulations 
with the various solvation terms. Fig. lb shows the corresponding plots for the 
susceptibility x(T) defined by 

X (T) =< n 2 H (T) >-< n H (T) > 2 . (16) 
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In Fig. la and lb the curves, representing the various simulations, fall into 
three groups. For the case where the protein-solvent interaction was approximated 
by a distance-dependent permittivity (DDE), both < uh > and x have a similar 
temperature dependence than is observed for poly-alanine in gas-phase simulations 
(GP). However, the transition temperature T c is shifted from T = 435 ± 20 K (gas- 
phase) to a higher value T = 495 ± 20. This temperature was determined from the 
maximum of the susceptibility x(T) in Fig. lb and is listed in table 1. To the same 
group belong the simulations in which the solvation energy was approximated by 
a solvent accessible surface term with either the OONS || or SCH [|l|] parameter 
set. In both cases susceptibility \ an d order parameter < riuiT) > show also a 
temperature dependence similar to the one of gas-phase simulations. Only now, 
the transition temperature T c is shifted to lower temperatures. The corresponding 
transition temperatures can be again determined from the positions of the maximum 
in x{T) and are also listed in table 1. The shift towards lower temperatures was one 
of the main results reported in Refs. [0, || for simulations with the OONS solvation 
energy, and our T c = 345 ± 20 K agrees well with their value T c = 340 K (no errors 
quoted) in Refs. 0, §• 

A somehow different behavior is observed in the simulation where the protein- 
water interaction was approximated by a solvent accessible surface term relying on 



the W92 [TT| parameter set. Here, the form of < rijj > indicates only partial helix 
formation and occurs only at much lower temperatures. The susceptibility x{T) in 
Fig. lb gives no indication for a helix-coil transition. For this reason no value of T c 
is listed for the W92 parameter set in table 1. Instead, we observe in Fig. 2 for this 
case at low temperatures even the appearance of residues whose backbone dihedral 
angles (ft, ip take values typical for a /3-sheet ( — 150 ± 30, 150 ± 30). 

Yet another behavior is observed in simulations where the solvation energy of 
Eq. is evaluated by means of the JRF parameter set. No formation of helices 
or sheets is observed in Figs. 1 and 2. Since no transition temperature can be 
determined, we do not list a value of T c for the JRF parameter set in table 1. 

The same grouping can be found in Fig. 3a-f where we display various energy 
terms as a function of temperature. In these figures we have shifted the solvation 
energies and the partial ECEPP/2 energies Eq, E lj , E hb and E tor of Eq. |5] by a 
constant term such that we have for all solvation models at T = 1000 K E so i = and 
Eq = Elj = Ehb = E tor = 0. Such a shift by an irrelevant constant allows a better 
comparison of the different simulations. The average total energy < E to t > which 
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is the sum of intramolecular potential energy Eeceep/2 and the solvation energy 
E so i, is displayed in Fig. 3a. We observe again that simulations with the parameters 
sets OONS and SCH, and such with distant dependent permittivity (DDE), have 
a similar temperature dependence as gas phase simulations (GP). On the other 
hand, in simulations relying on the W92 parameter set, the energy varies less with 
temperature and is at low temperatures considerably higher than in the simulations 
with other solvation energy terms. Finally, the energy in simulations with the JRF 
parameter set is an almost linear function of temperature and is especially at high 
temperatures much lower than the energies found in gas phase simulations. 

The dissimilar behavior of energy for simulations with different solvation terms 
is even more obvious in Fig. 3b where the average intramolecular energy Eecepp/2 
is drawn. While this energy term decreases between 1000 K and 150 K by ~ 50 
Kcal/mol (with most of that change, ~ 30 Kcal/mol, happing around the respective 
transition temperature T c ) in gas-phase simulations (GP) and in simulations with 
OONS, SCH and DDE solvation terms, it changes in the same temperature interval 
only by w 20 Kcal/mol in simulations utilizing the JRF or W92 parameter sets. 
Since for these two parameter sets also no or only little helix formation was observed 
it seems likely that the formation of helices is related to the large gain in potential 
energy observed for GP, OONS, SCH and DDE simulations. This gain in potential 
energy is in part due to the formation of hydrogen bonds between a residue and 
the fourth following one in the polypeptide chain which stabilize an a-helix. Fig. 3c 
displays the average hydrogen-bonding energy < E HB > of Eq. || as a function of 
temperature and one can clearly see the gain in energy for the GP, DDE, OONS and 
SCH simulations at the respective helix-coil transition temperatures of table 1. No 
such gain is observed in W92 and JRF simulations where also no helix formation was 
found. A similar gain in energy with helix formations in gas-phase and simulations 
with DDE, OONS and SCH solvent representations is also observed for the average 
Lennard- Jones energy < E LJ > and the electrostatic energy < E c > displayed in 
Fig. 3d and 3e, respectively. Note also in Fig. 3e the large gain in Ec for DDE at 
the helix-coil transition temperature which additionally stabilizes the a-helix in this 
model. 

A complementary picture is found in Fig. 3f where the solvation energy E so i 
is shown as a function of temperature. The observed helix formation in gas phase 
simulations and such with OONS, SCH and DDE solvent representations is correlated 
with an increases of the solvation energies by ~ 5 Kcal/mol. On the other hand, 
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in simulations with the W92 and JRF parameter sets, for which no helix-formation 
was observed in Fig. 1, E aoi decreases with temperature. This decrease is only rj 5 
kcal/mol for W92, but it is much larger (of order 30 kcal/mol) in the case of JRF 
where the solvation energy is the dominant term. 

The effects of the dominant solvation term in simulations with the JRF parameter 
set can also be seen in Fig. 4. In this figure the average radius of gyration, a measure 
for the compactness of configurations, is shown as a function temperature. One can 
see that this quantity changes little with temperature for the JRF data. However, 
its value is over the whole temperature range considerably smaller than observed in 
the other simulations. This indicates that the JRF solvation term favors already at 
high temperatures compact configurations, and that the pressure towards compact 
structure is such that the more elongated helices cannot be formed. Note however, 
that the tendency towards compact configurations does not lead to a lower Lennard- 
Jones energy E LJ as one would expect. Fig. 3d indicates that < E LJ > is at low 
temperatures in JRF simulations even larger than in GP, DDE, OONS and SCG 
simulations where helix-formation was observed. The tendency towards compact 
structures in JRF simulations may be due to the fact that JRF parameter set was 
developed from minimum energy (i.e. compact) conformations of peptides (the 
low-energy conformations of 13 tetrapeptides derived by NMR studies pD|), and 
therefore this parameter set may have an intrinsic bias towards compact structures. 

On the other hand, the W92 parameter set was developed from measurements of 
free energies of amino acid side chain analogs from vapor to water |2(J . The param- 
eters for this set are negative for all atoms except carbon meaning that the nitrogen, 
oxygen and sulfur atoms are considered hydrophilic, i.e. favoring solvent exposure. 
This explains not only the small solvation energies observed for this parameter set 
in Fig. 3f, but also why in Fig. 4 the radius of gyration is consistently larger for 
this parameter set than for the others indicating that extended configurations are 
favored with this parameter set. This bias towards extended structures limits again 
the formation of a-helices. 

While the OONS parameter set was derived from experimental free energies of 
gas-to-water transfer of small aliphatic and aromatic molecules, the SCH is not 
directly based on experimental free energy data. Instead, it was developed as an 
optimized parameter set to complement the CHARMM force field p5|] . In both 
parameter sets the hydrophobic character of the carbon atoms is increased and the 
hydrophilic character of uncharge oxygen and nitrogen atoms decreased resulting 
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into the large solvation energies of these two parameter sets (when compared with 
the one of the W92 parameter set) that one observes in Fig. 3f. The OONS and SCH 
solvation energies again favor extended structures (the radius of gyration has larger 
values than found in gas-phase simulations), however, the interplay of solvation 
energies and intramolecular ECEPP/2 energy is such that the radius of gyrations 
(and consequently the compactness) of polyalanine configurations as a function of 
temperature shows a similar behavior as the gas-phase simulation. The same is true 
for the DDE simulation where the protein-solvent interaction was approximated by 
a distance- dependent permittivity. 

Our results demonstrate that the helix formation is due to the gain in poten- 
tial (intramolecular) energy while (with the exception of the JRF parameter set) 
the solvent-accessible surface terms favor extended peptide configurations. Table 2 
summarizes the differences in total energy AE tot , solvation energy AE soh poten- 
tial energy AE ECEE p/2, and the partial energies AE C) AE LJ) AE HB and AE tor 
between complete helical configurations (all residues with exception of the terminal 
ones are part of an a-helix) and coil configurations at temperature T = 300 K for 
gas-phase, DDE OONS and SCH simulations. Note, that the intramolecular energy 
differences AE EC epp/2 of gas-phase, OONS and SCH simulations have within their 
error bars the same values. For simulations with the W92 parameter set the longest 
found helix consists of 6 consecutive residues. Hence, we measured for this case only 
the energy difference between configurations with at least three consecutive helical 
residues (i.e. one turn of an a-helix) and coil configurations. This modified defini- 
tion of the energy differences is also the reason for the smaller value of AE EC epp/2 
listed for W92 in table 2. We do not list energy differences for the JRF parameter 
set since no helices were found in simulations utilizing this parameter set. 

Note that in simulations with distant dependent permittivity (DDE) helices are 
energetically more favored than in the gas-phase simulations. This is due to the 
increased contribution from the Coulomb term E c as one can also see in Fig. 3e. 
The larger energy gap between helical and coil conformations (when compared with 
gas-phase simulations) explains why the transition temperature is higher in DDE 
simulations than in gas-phase simulations. 

For the OONS and the SCH parameter set the solvation energy difference AE so[ is 
positive (indicating that coil structures are energetically favored), but its magnitude 
is only approximately half that of the potential energy difference AE EC epp/2- Hence, 
there is still an overall energetic gain connected with helix formation. However, 
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in both cases the total energy difference between helical and coil configurations is 
reduced by the solvation energy when compared with the gas-phase simulation. This 
reduction of the energy gap leads to the lower transition temperatures observed in 
OONS and SCH simulations. 

On the other hand, for the W92 parameter set we find that AEeccp/2 and 
AE so i are of same magnitude so that helical configurations are not or only weakly 
energetically favored. This is consistent with our results in Fig. la and lb where 
we find at T = 280 K a high average helicity in OONS and SCH simulations but 
only a small value of < n# > and no indications for a helix-coil transition in W92 
simulations. An evaluation of energy differences was not possible for simulations 
with the JRF parameter set since no helices were found. 

The above results indicate that the existence and characteristics of the helix-coil 
transition in polyalanine depend strong on the details of the solvent representation. 
In order to evaluate the validity of the different solvent models one has to compare 
the numerical results with experimental data. For this purpose we have calculated 



the helix propagation parameter s which was also determined by experiments [28 



p9[ . According to the Zimm-Bragg model |27] the average number of helical residues 
< n > and the average length < i > of a helical segment are given for large number 
of residues N by 

1 1-5 



N 2 2.^(1 - s) 2 + Asa 

2s 



■ 



<i> = 1+ 



1- s + J(l- s) 2 + Asa 



where s is the helix propagation parameter and a the nucleation parameter of the 
Zimm-Bragg model. From these equations with the values of < n >/N and < I > 
calculated from the multicanonical production runs, we have calculated s at tem- 
perature T = 280 K for gas-phase and the different solvation models. Our values are 
summarized in table 3 which also lists our a values. Our results for gas-phase, DDE 



and OONS simulations are in agreement with the experimental results of Ref. [59 
where they list values of s(Ala) between 1.5 and 2.19. On the other hand, the s 
value obtained in the SCH simulation agrees well with the one obtained by the host- 



guest technique of Ref. [ 2S| . However, the s values which were obtained in W92 or 
JRF simulations do not agree with either of the experimental data. Hence, we con- 
clude that the W92 and JRF parameter sets are not appropriate solvation models 
in simulations of polyalanine. Otherwise, the variation in the experimental data is 
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too large to give indications whether one of the remaining solvent representations 
(DDE, OONS, SCH or even no solvent at all (GP)) is preferable over the others. 

4 Conclusions 

We have performed multicanonical simulations of polyalanine. The intramolecular 
forces were modeled by the ECEPP/2 force field and various approximations for the 
solvation energy were studied. We observed that whether a helix-coil transition is 
observed for poly-alanine, and at what temperature, depends strongly on the chosen 
approximation for the protein-solvent interaction. Our results demonstrate both the 
importance (and need) of including solvation terms into protein simulations and the 
difficulties in chosing an adequate representation of the protein-water interactions. 
Especially when using the solvent-accessible surface approach, it seems necessary to 
carefully choose a parameter set that is suitable for the problem under consideration. 
Use of a specific parameter set without further justification could otherwise generate 
miss-leading results. 

Acknowledgement 

U. Hansmann gratefully acknowledges support by a research grant from the National 
Science Foundation (CHE-9981874). This article was written in part while U.H. was 
visitor at the Department of Physics at University of Central Florida. He thanks 
Alfons Schulte, Weili Luo, Aniket Bhattacharya and Brian Tonner for their kind 
hospitality during his stay in Orlando. 

References 

[1] D. Poland and H.A. Scheraga, Theory of Helix- Coil Transitions in Biopolymers 
(Academic Press, New York, 1970). 

[2] Kemp, J. P., and Chen,Z.Y. 1998. Formation of Helical States in Wormlike Poly- 
mer Chains. Phys. Rev. Lett. 81:3880-3883. 

[3] Hansmann, U.H.E., and Okamoto Y. 1999. Finite-size scaling of helix-coil tran- 
sitions in poly-alanine studied by multicanonical simulations. J. Chem. Phys. 
110:1267-1276; 111: 1339(E). 



11 



[4] Alves, N.A., and Hansmann, U.H.E. 2000. Partition Function Zeros and Finite 
Size Scaling of Helix-Coil Transitions in a Polypeptide. Phys. Rev. Lett. 84: 
1836-1839. 

[5] Alves, N.A., and Hansmann, U.H.E. 2001. Yang-Lee zeros and the helix-coil 
transition in a continuum model of polyalanine. Physica A 292: 509-518. 

[6] Hudgins,R.R., Ratner, M.A., and Jarrold, M.F. 1998. Design of Helices that 
are stable in vacuo. J. Am. Chem. Soc. 120: 12974-12975. 

[7] Mitsutake, A., and Okamoto, Y. 1999. a-Helix propensities of homo-oligomers 
in aqueous solution my multicanonical algorithm. Chem. Phys. Let. 309: 95- 
100. 

[8] Mitsutake, A., and Okamoto, Y. 2000. Helix-coil transitions of amino-acid 
homo-oligomers in aqueous solution studied by multicanonical simulations. J. 
Chem. Phys. 112: 10638-10647. 

[9] Ooi,T., Obatake,M., Nemethy, G., and Scheraga, H.A. 1987. Accessible surface 
areas as a measure of the thermodynamic parameters of hydration of peptides. 
Proc. Natl. Acad. Set. USA 8:3086-3090. 

[10] Vila, J., Williams, R.L, Vasquez, M., and Scheraga, H.A. 1991. Empirical solva- 
tion models can be used to differentiate native from near-native conformations 
of bovine pancreatic trypsin inhibitor. Proteins Struct Fund Genet 10: 199-218. 

[11] Wesson, M., and Eisenberg, D. 1992. Atomic solvation parameters applied to 
molecular dynamics ofproteins in solution. Protein Science 1:227-235. 

[12] Schiffer, C.A., Caldwell, J.W., Kollman, PA., and Stroud, R.M. 1993. Protein 
structure prediction with a combined solvation free energy molecular mechanics 
force field. Mol. Simul. 10: 121-149. 

[13] Sippl, M.J., Nemethy, C, and Scheraga, H.A. 1984. Intermolecular potentials 
from crystal data. 6. Determination of empirical potentials for O-H- • -0=C 
hydrogen bonds from packing configuration. J. Phys. Chem. 88: 6231-6233; 
and references therein. 

[14] Eisenmenger, F., Hansmann, U.H.E., Hayryan,Sh., and Hu, C.K. 2001. [SMMP] 
A Modern Package for Simulation of Proteins. Comp. Phys. Comm. 138: 192- 
212. 



12 



[15] Hingerty, B., Richie, R.H., Ferrel, T.L., and Turner, J. 1985. Dielectric Effects 
in Biopolymers-The Theory of Ionic Saturation Revisited. Biopolymers 24: 
427-439. 

[16] LaveryR., Sklenar, H., Zakrzewska, K., and Pullman, B. 1986. The Flexibility 
of the Nucleic Acids: (II) The Calculation of Internal Energy and Applications 
to Mononucleotide Repeat DNA. J. Biomol. Struct. & Dynamics 3: 989-1014. 

[17] Lee, B., and Richards, F.M. 1971 Interpretation of Protein Structures: Estima- 
tion of Static Accessibility. J. Mol. Biol. 55(3): 379-400. 

[18] Eisenberg, D., and McLachlan, A.D. 1986. Solvation energy in protein folding 
and binding. Nature 319: 199-203. 

[19] Berg, B.A., and Neuhaus, T. 1991. Multicanonical algorithms for first order 
phase transitions. Phys. Lett. B 267:249-253. 

[20] Hansmann, U.H.E., and Okamoto, Y. 1998. The Generalized-Ensemble Ap- 
proach for Protein Folding Simulations, in: Stauffer, D. (ed.) " Annual Reviews 
in Computational Physics VP (Singapore: World Scientific), 129-157. 

[21] Hansmann, U.H.E., and Okamoto, Y. 1993. Prediction of peptide conformation 
by multicanonical algorithm: a new approach to the multiple-minima problem. 
J. Comp. Chem. 14: 1333-1338. 

[22] Okamoto, Y., and Hansmann, U.H.E. 1995. Thermodynamics of Helix - Coil 
Transitions Studied by Multicanonical Algorithms. J. Phys. Chem. 99: 11276- 
11287. 

[23] Ferrenberg, A.M., and Swendsen, R.H. 1988. New Monte Carlo technique for 
studying phase transitions. Phys. Rev. Lett. 61:2635-2638. 

[24] Berg, B.A. 1996 Multicanonical Recursions. J. Stat. Phys. 82: 323-342. 

[25] Brooks, B.R., Bruccolieri, R.E., Olafson, B.D., States, D.J., Swaminathan, 
S., and Karplus, M. 1983. CHARMM: a program for macromolecular energy, 
minimization, and dynamics calculations. J. Comp. Chem. 4: 187- 217. 

[26] Wolfenden, R., Andersson, L., Cullis, P.M., and Southgate, C.C.B. 1981. Affini- 
ties of Amino-Acid Side Chains for Solvent Water. Biochemistry 20: 849-855. 



13 



[27] Zimm B.H., and Bragg, J.K. 1959. Theory of the Phase Transition Between 
Helix and Random Coil in Polypeptide Chains. J. Chem. Phys. 31: 526 -535. 

[28] Wojcik, J., Altmann, K.H., and Scheraga, H.A. 1990. Helix-Coil Stability Con- 
stants for the Naturally-Occuring Amino-Acids in Water. 24. Half-Cystine 
Parameters from Random Poly(Hydroxibutylglutamnie CO-S-Methylthiol- 
L_cystine) Biopolymers 30, 121-134. 

[29] Chakrabartty, A.; R.L. Baldwin, R.L. 1993. In: Protein Folding: In Vivo and 
In Vitro; Cleland J.; King, J. eds.; ACS Press: Washington, D.C., pp. 166-177. 



14 



Tables: 



Table 1: Transition temperatures for the helix-coil transition in ALA 10 as obtained 
from gas-phase simulations and simulations with various solvent representations. All 
results rely on multicanonical simulations of 1,000,000 Monte Carlo sweeps each. 



Model 


T 


GP 


435(20) 


DDE 


495(20) 


OONS 


345(15) 


SCH 


285(25) 


W92 




JRF 





Table 2: Energy differences between helical and configurations (see text) at T = 
280 K as measured in gas-phase simulations and simulations with various solvent 
representations. All results rely on multicanonical simulations of 1,000,000 Monte 
Carlo sweeps of ALA 10 for each case. 



Model 


A£ tot 


AE sol 


^E E cEPP/2 


AE C 


A£ LJ 


AE HB 


AE tor 


GP 


-16.9(1) 




-16.9(1) 


0.4(3) 


-12.1(1) 


-4.3(3) 


-0.8(1) 


DDE 


-17.9(6) 




-17.9(6) 


-3.6(2) 


-10.1(4) 


-3.9(2) 


-0.3(1) 


OONS 


-11.3(9) 


4.1(3) 


-15.4(6) 


-0.2(1) 


-10.7(4) 


-4.1(1) 


-0.4(1) 


SCH 


-7.1(5) 


8.7(1) 


-15.8(5) 


0.7(3) 


-11.2(2) 


-4.6(3) 


-0.7(1) 


W92 


-0.7(7) 


5.6(6) 


-6.3(1.1) 


0.8(1) 


-5.8(9) 


-1.0(2) 


-0.3(1) 


JRF 
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Table 3: Helix propagation parameter s and nucleation parameter a at T = 280 K 
for Ala 10 as measured in gas-phase simulations and simulations with various solvent 
representations. All results rely on multicanonical simulations of 1,000,000 Monte 
Carlo sweeps for each case. 



Model 


s 


a 


GP 


1.67(9) 


0.15(1) 


DDE 


1.78(12) 


0.15(1) 


OONS 


1.31(15) 


0.13(1) 


SCH 


1.02(15) 


0.11(2) 


W92 


« 


> 1 


JRF 


« 


> 1 
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Figure Captions: 



Fig. 1 Temperature dependence of (a) the average number < n H > of helical 
residues and (b) the susceptibility x(T) for ALA 10 as calculated from a gas- 
phase simulation and from simulations with various solvation energy terms. 
All results rely on multicanonical simulations of 1,000,000 Monte Carlo sweeps 
each. 

Fig. 2 Temperature dependence of the average number < % > of residues whose 
backbone dihedral angles 0, ip take values as typically found in /3-sheets. Re- 
sults from a gas-phase simulation and such with various solvation terms are 
presented for ALAio- All data rely on multicanonical simulations of 1,000,000 
Monte Carlo sweeps. 

Fig. 3 Temperature dependence of (a) the total energy < E tota i = E E ceep/2 + 
E S oi >, (b) the intramolecular energy < E EC epp/2 >, (c) the hydrogen- 
bonding energy < E HB >, (d) Lennard- Jones energy < E LJ >, (e) Coulomb 
energy < E c >, and (f) the solvation energy < E so i > as calculated from 
a gas-phase simulation and from simulations with various solvation energy 
terms. All results rely on multicanonical simulations of ALAi with 1,000,000 
Monte Carlo sweeps for each case. 

Fig. 4 Temperature dependence of the average radius-of-gyration < R gy > as mea- 
sured in gas-phase simulations and simulations with various solvent represen- 
tations. All data rely on multicanonical simulations of 1,000,000 Monte Carlo 
sweeps. 
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